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Abstract 


The vertical structure equation is a singular Sturm- Liouville problem 
whose eigenfunctions describe the vertical dependence of the normal modes of 
the primitive equations linearized about a given thermal profile. The 
eigenvalues give the equivalent depths of the modes. We study the spectrum 
of the vertical structure equation and the appropriateness of various upper 
boundary conditions, both for arbitrary thermal profiles. Our results 
depend critically upon whether or not the thermal profile is such that the 
basic state atmosphere is bounded. 

In the case of a bounded atmosphere we show that the spectrum is 
always totally discrete, regardless of details of the thermal profile. For 
the barotropic equivalent depth, which corresponds to the lowest eigenvalue, 
we obtain upper and lower bounds which depend only on the surface 
temperature and the atmosphere height. All eigenfunctions are bounded, but 
always have unbounded first derivatives. We prove that the commonly invoked 
upper boundary condition that vertical velocity must vanish as pressure 
tends to zero, as well as a number of alternative conditions, is well posed. 

For unbounded atmospheres, on the other hand, we show that typically 
there is a continuous spectrum, the boundary condition of vanishing vertical 
velocity is not well posed, and the eigenfunctions, if any, are unbounded. 
We point out, however, that owing to the traditional shallowness 
approximations of the primitive equations, the vertical structure equation 
has no meaning for unbounded atmospheres. This leads to the conclusion that 
the vertical structure equation always has a totally discrete spectrum under 
the assumptions implicit in the primitive equations. 
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1. Introduction and Main Results 

Atmospheric normal modes play a central role in the theory of tidal 
motions of the atmosphere (Siebert, 1961; Chapman and Lindzen, 1970) and in 
the theory of planetary- scale traveling waves (Salby, 1984). Normal modes 
are also important in numerical weather prediction, where they are used to 
perform objective analysis (Flattery, 1970), to provide basis functions for 
spectral models (Kasahara, 1977), to initialize models (Daley, 1981), to 
represent global data (Kasahara and Puri, 1981), and to construct test 
solutions for models (Dee and da Silva, 1986). Recently normal modes have 
also been used as the basis of a stability theory for four -dimensional data 
assimilation methods (Cohn and Dee, 1986). 

The vertical structure equation (VSE) is a Sturm- Liouville ordinary 
differential equation whose eigensolutions describe the vertical structure 
of the normal modes of the linearized primitive equations. Since the VSE is 
singular at the top boundary, where pressure vanishes, two fundamental 
questions are: 

(A) What is the nature of the spectrum of the VSE, i.e. , when is 
the spectrum totally discrete and when is there a continuous spectrum? 

(B) What boundary condition is appropriate at the top boundary? 

The basic state thermal profile appears in a coefficient of the VSE, and 
various authors have considered these two questions for particular choices 
of the thermal profile (Taylor, 1936; Dikii, 1965; Wiin-Nielsen, 1971; 
Staniforth et al . . 1985). Our objective in this paper is to answer these 
questions for arbitrary profiles. 

Before summarizing results, we review the background and derivation of 
the VSE following Daley (1981) and Staniforth et al . (1985). In the usual 
rotating pressure coordinate system, the adiabatic, hydrostatic primitive 

equations linearized about a state of rest can be written as 
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Here u, v and w are the perturbation eastward, northward and vertical 
velocity components in pressure coordinates, and $ is the perturbation 
geopotential. The independent variables are time t, pressure p, longitude A 
and latitude The constants a and Q are the earth's radius and rotation 

rate. The thermal structure of the basic state enters through the static 
stability function a. 
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where Tq — Tq(p) is the basic state temperature profile and 8 ^ 


T 0 p * is 


the corresponding potential temperature profile; R is the gas constant for 
dry air and k - R/c^, c^ being the specific heat at constant pressure. 

We take the domain of (1.1) to be global, and with pressure ranging 
from a constant surface value p - p^ up to p - 0. The lower boundary 
condition for (1.1) corresponding to no mass flux through the earth's 


surface is 
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We do not yet specify an upper boundary condition. Note, however, that the 
time rate of change of total energy is given by 
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where dA is the horizontal area element, the triple integral is over the 
whole domain, and the double integral is over a spherical shell. Thus a 
boundary condition 


-* 0 as p ■* 0 


( 1 . 5 ) 


would guarantee energy conservation, while the commonly invoked condition 


w -» 0 as p •* 0 


( 1 . 6 ) 


does not necessarily conserve energy, because one does not know a priori 
whether $(t,A,^,p) is bounded for all time. 

Seeking separable solutions 
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one obtains the vertical structure equation and lower boundary condition 


(LBC) 
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along with the linear global shallow-water system 
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whose solutions are the Hough harmonics (Hough, 1898; Longuet-Higgins, 1968; 
Holl, 1970; Kasahara, 1976). Equations (1.8) and (1.10) are coupled through 
the separation constant n which appears as an eigenvalue in the VSE (1.8). 
The corresponding equivalent depth H is H - 1/fig, g being the acceleration 
uue to gravity. The vertical velocity w is obtained diagnostically as 
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The candidate boundary conditions (1.5) and (1.6) imply, respectively, that 


; X g -* 0 as p -♦ 0 , (1.12) 

or 

0 

a dp 


as p -► 0 . 


(1.13) 
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Staniforth et al . (1985) studied the VSE for the two-parameter family 
of static stability functions 


a(p) - Ap - * , (1.14) 

where A > 0 and a are constants , which includes many of the static stability 
functions considered by previous authors. They showed that for the linear 
system (1.1) to be valid in the sense of the neglected nonlinear terms being 
smaller than the retained linear terms, it is necessary that a > 2— k. On 
the other hand, for Tg(p) to be finite as p -* 0, they showed that a :S 2. 
Hence they assumed that 


2-k s a :S 2 . (1.15a,b) 

In case a < 2, they were able to solve the VSE explicitly in terms of Bessel 
functions. These solutions were used to show that the spectrum is totally 
discrete, i.e., there is an infinite sequence of eigenvalues 

^0 < M 1 < ^2 < ' ' * * * ( 1 - 16 ) 

for which there is a corresponding sequence of eigenfunctions x n satisfying 
the VSE (1.8), the LBC (1.9), and the upper boundary condition (1.13). For 
a - 2, which includes the much-studied isothermal case, they showed that 
there is at most one eigenvalue with corresponding eigenfunction satisfying 
(1.8,1.9,1.13), while the rest of the spectrum is continuous. 

The VSE is a s insular S turra - Li ouvi 1 1 e equation whenever the 
coefficient o(p) ■* « as p -* 0; in particular it is singular for the family 
of profiles (1.14,1.15). The answer to questions (A) and (B) , as well as 
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other qualitative aspects of solutions, depend solely on the rate at which 
a(p) •* <*. More precisely, they depend only on the behavior of o(p) , and 
hence on that of T^(p) , in an arbitrarily small neighborhood of the singular 
point p - 0. For example, we will see that by modifying the isothermal 
profile near p - 0, the continuous spectrum becomes totally discrete. 

In analyzing the VSE, one must ensure that any assumptions on the 
asymptotic behavior of Tq(p) are consistent with the assumptions implicit in 
the governing nonlinear primitive equations. Recall that the primitive 
equations are based on certain "shallowness" or "traditional" approximations 
(Phillips, 1966, 1973), in which gravity g (the magnitude of the resultant 
of the centrifugal and Newtonian gravitational forces) is taken to be 
constant and in which the metric factors are replaced by their values near 
the surface of the earth. These approximations break down for unbounded 
atmospheres. Therefore we can be assured that an analysis of the linearized 
system (1.1), and therefore of the VSE, is meaningful only when the basic 
state atmosphere is bounded. 

Boundedness of the basic state atmosphere can in fact be expressed 
as a condition on Tq(p) near p - 0, as follows. Since the basic state is 
hydrostatic, the basic state geopotential $q is given by 

?S RT n (s) 

* 0 (p) - W + J “V" ds ; (l - l7) 

P 

cf. Staniforth et al . (1985, eq. 2.17). Hence the geopotential extent of 
the domain is 

I S RT ( s ) 

^ ds , (1.18) 
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and the basic state atmosphere is bounded if and only if the limits in 
(1.18) exist and are finite, i.e. 


?S RT n (s) 

$ E “ I g ds < * . (1.19) 

0 

In particular, if Tg(p) has a limit as p -*• 0 then a necessary condition for 
the basic state atmosphere to be bounded is that in fact 

T 0 (p) ^ 0 as p - 0 , (1.20) 

for otherwise the integral in (1.19) would diverge at least logarithmically. 
In a recent PhD thesis (Sudarshan, 1985) , the plausibility of condition 

(1.20) has also been argued for atmospheres in conductive equilibrium. 
Conditions (1.19) and (1.20) will be fundamental in our analysis. We 

will show in particular that (1.20) guarantees a totally discrete spectrum 
and that under the stronger condition (1.19), but not necessarily under 

(1.20) , the boundary conditions (1.12) or (1.13), amongst others, are indeed 
appropriate. 

One might argue that neither of conditions (1.19) and (1.20) is 
physically reasonable. For example, the temperature of the thermospheric 
plasma varies from ^bout 600K to as much as about 2000K, depending on the 
amount of solar activity (Gill, 1982, pp. 48-49). Nevertheless, our point 
of view is that in writing down the vertical structure equation one is 
addressing the primitive - equation model of the atmosphere, not the 
atmosphere itself. The shallowness assumptions of the primitive equations 
then make discussion and analysis of unbounded atmospheres irrelevant. In 
this paper we do give some results for unbounded atmospheres, partly in 
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Theorem 1 and partly by way of examples in Sec. 4, but only for the sake of 
contrast with our results for bounded atmospheres. 


The following basic assumptions are made throughout our analysis. 
Static stability assumptions . For 0 < p < p g we asstime that 


and 


It follows from (1.2) that 


dp 


< 0 


T q > ° . 


( 1 . 21 ) 

( 1 . 22 ) 


0 < cr (p) < « for 0 < p S p s , (1.23) 


so that the VSE is regular on any closed interval not containing p «• 0. 

Continuity assumptions . We assume that a and da/dp are continuous for 
0 < p ^ Pg- Together with (1.23), this guarantees that all of the classical 
theory of singular Sturm- Liouville problems applies. 

Limit assumptions . We make the technical assumption that the basic 

dT 0 

state temperature profile Tq(p) is such that both Tq and have limits 

as p -* 0. Thus we exclude certain pathological profiles that oscillate 
infinitely often near p - 0. We also asstime that lim T n (p) is finite. 

p-0 u 

Entropy assumptions . We assume that 


Also we assume that 


lim inf ( 

p-*0 



d0 ( 

dp 


) > 0 . 


(1.24) 


1/4 < k < 1/2 ; 


(1.25) 
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for dry air, k - 2/7 for temperatures in the range 250-400K (Batchelor, 
1967, p. 45). In Sec. 2 it will be shown that (1.24) implies that the 
singularity of the VSE is of limit-point type . There it will also be seen 

2+k 

that (1.24) holds whenever 8^ •* <*> and p dfl^/dp has a limit, possibly 
infinite, as p -*■ 0. In particular, (1.24) is satisfied for static stability 
functions of power-law form (1.14) when the linearization validity criterion 
(1.15a) is met. 

Our first main result, proven in Sec. 2, is the following theorem. 


Theorem 1 . Consider the VSE . (1.8), LtfC (1.9) and upper boundary 
condition 


and let 



(1.26) 
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00 


lim T 0 (p) . 
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If T - 0 then the spectrum is totally discrete, i.e., there is a countably 
infinite set of eigenvalues [n n ) , 

p 0 < M 1 < M 2 < "• “ 40 ’ (1-28) 

and corresponding eigenfunctions {X^} for which (1.8, 1.9, 1.26) are 
satisfied. If T^ * 0 then the portion of the spectrum for 

a < u =■ - — ~ (1.29) 

M M c 4/cRT 


is discrete. 
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The basic state atmosphere must be bounded for consistency with the 
shallowness approximations in the nonlinear primitive equations. This, in 
turn, implies that T =0. Theorem 1 therefore allows us to conclude that 
under the assumptions implicit in the primitive equations, the VSE always 
has a totally discrete spectrum. 

In general the square- integrability upper boundary condition (1.26) is 
the only one which is well-posed in the sense of always excluding one of the 
two linearly independent solutions of the VSE for every p. The following 
theorem shows that for bounded basic state atmospheres, however, there are a 
number of equivalent upper boundary conditions. This theorem is proven in 
Section 3. 


Theorem 2 . Consider the VSE (1.8) subject to the LBC (1.9). If the 
basic state atmosphere is bounded, i.e. if (1.19) holds, then all of the 
following are equivalent upper boundary conditions : 


(i) 



X 2 (p) dp < « 


(ii) 

(iii) 


lim |X(p) | < «* 

p->0 


lim 
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n a dp 
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(1.30a) 
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(1.30c) 


(1 . 30d) 


lim p^X - 0 

p-0 


(v) 


if 0 < 5 < 1/2 . 


(1 . 30e) 


The equivalence of conditions (i) and (ii) implies that the eigen- 
functions for a bounded atmosphere are bounded. It also follows from the 
theorem that for bounded atmospheres the boundary conditions (1.12, 1.13) 
implied by -» 0 (1.5) and w -*• 0 (1.6), respectively, are in fact 
equivalent and well posed. Any of conditions (ii)-(v) could be useful in 
numerical calculations . 

In Section 3 we also prove the following theorem, which summarizes our 
main results for bounded atmospheres . 


Theorem 3 . Consider the VSE (1.8) for a bounded atmosphere, subject 
to the LBC (1.9) and any one of the equivalent upper boundary conditions of 
Theorem 2 . Then 

(i) The spectrum is totally discrete. 

(ii) The eigenvalues are all positive and satisfy 



(1.31a) 


*« e - rt s < " n 


for n > 1 , 


(1.31b) 


where is the geopotential extent defined in (1.18) and Tg is the surface 

temperature. In terms of the equivalent depths - (g/i^). \ these 
inequalities are 

RT S < gH Q < k$ e (1.32a) 


gH < * (/c$i -RT ) 
b n n ' E S 


n > 1 . 


(1.32b) 


14 


(iii) The eigenfunctions X^(p) form a complete orthogonal basis for 
the space of functions which are square- integrable over 0 < p < p_. 
Moreover, X n has precisely n zeros on 0 < p < p g . 

(iv) While each eigenfunction X^ remains bounded as p -*• 0, the 
derivative is unbounded. In fact there is a constant c - c(n) > 0 such that 


dX 
n 

dp 


> cp 


- 1/2 


(1.33) 


for all p sufficiently small. 

Statement (i) of this theorem follows directly from Theorem 1. The 
bounds in statement (ii) are independent of the details of the thermal 
profile: they depend only on the surface temperature T g and geopotential 
extent in addition to the constants k , R, and g. For shallow 
atmospheres, /c$g is not much larger than RT g . The bound (1.32a) therefore 
implies that the barotropic equivalent depth Hq is not sensitive to details 
of the thermal profile. 

Statement (iii) justifies the use of the vertical structure functions 
X n as basis functions for spectral expansions. Statement (iv) implies that 
there is a "boundary layer" near p - 0, where each vertical structure 
function changes rapidly. Numerical methods for both the VSE and the 
primitive equations themselves should be designed to take this into account. 

If one assumes an asymptotic form for T Q (p) near p - 0 then the 
asymptotic form of the eigenfunctions can be obtained. This is done in 
Sec. 4. As a first example we consider 

00 

T (p) = 2 r p €J 
j-0 J 


near p - 0 , 


(1.34) 
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where 0 < e < k, r Q - T^ - 0, and r \ > Suc ^ an atmosphere is bounded, 
and the equivalence of the boundary conditions in Theorem 2 is verified by 
displaying the asymptotic behavior of the eigenfunctions. 

Second, we consider unbounded atmospheres for which 

ao 

T 0 (P) - + 2 r pj near p - 0 , (1.35) 

with T > 0, which includes all atmospheres with isothermal tops. In this 
case we extend the result (1.29) of Theorem 1 by using the asymptotic 
behavior of the eigenfunctions to show that the spectrum is in fact 
continuous for n S n . Hence there is only a finite number of eigenvalues 
in this case. 

For further contrast we consider two families of atmospheres for which 
Tq(p) is specified throughout the domain but neither (1.34) nor (1.35) is 
satisfied. For both families Tg(p) ■* 0 as p ■* 0, but so slowly that the 
atmosphere is unbounded. We exhibit a closed-form solution which is square- 
integrable but unbounded, and another which is not square -integable yet has 
^ ^ 0 as p -*• 0. Hence the hypothesis of Theorem 2 that the basic state 

atmosphere be bounded cannot be weakened. 

We conclude with a brief discussion of our results in Sec. 5. 
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2 . Spectral theory 

The object of this section is to prove Theorem 1, which shall be done 
by application of the classical theory of singular Sturm- Liouville problems. 
Of fundamental importance is that the VSE is always of limit-point type, 
which we will first elaborate upon. 

From general existence theory it is known that the VSE (1.8), with no 
boundary conditions imposed, has precisely two linearly independent 
solutions for any number /i. A Sturm- Liouville differential equation is said 
to be of limit-point type if for some choice of fi, say p - Hq, one solution 
is not square- integrable. The classical theory due to Weyl (cf. Coddington 
and Levinson, 1955, Ch. 9, Theorem 2.1) shows that in fact such an equation 
has at most one square - integrable solution for anv number p. Thus an 
appropriate boundary condition for equations of limit-point type is square- 
integrability , since this condition guarantees uniqueness. The spectrum is 
then determined by imposing a boundary condition at the regular point, in 
our case the LBC (1.9). 

By choosing - 0, to demonstrate that the VSE is of limit-point type 
it suffices to show that the homogeneous equation 

d£ d^> “ 0 > o < p < p s , (2.1) 

has a solution which is not square- integrable. One solution is X - constant, 
which of course is square - integrable . The other solution is X - VKp) » 



P 


a(s) 


ds . 


( 2 . 2 ) 


We will show that is not square - integrable , i.e.. 
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r s 2 

i> (s) ds -* ® as p -*• 0 , 


(2.3) 


by proving the following lemma. 


Lemma 2.1 . There exists a constant > 0 such that 


a(p) > c^p 3/ ^ for 0 < p < p s . 


(2.4) 


Proof : The static stability assumption (1.21) implies that 


2+k d ^ 

-P dp 


> 0 for 0 < p S p. , 


(2.5) 


and by virtue of the entropy assumption (1.24) there is a constant S > 0 
such that in fact 


2 + * «0 „ . 

■P dT > ‘ 


for 0 < p ^ p< 


( 2 . 6 ) 


Substituting (2.6) into the definition (1.2b) of a, we then have 


■ o (p) > SRp 3/2 for 0 < p 3 p s . (2.7) 

Lemma 2.2 . The VSE is of limit-point type, so that an appropriate 
upper boundary condition is 

Ps 

| X 2 (s) ds < « . (2.8) 

0 



Proof : Inserting the bound (2.4) into the definition (2.2) of -0(p) , 

one obtains by integration 


*(p) > 2c^(p 1/2 - p s 1/2 ) . (2.9) 

Therefore the integral in (2.3) diverges at least logarithmically, and the 
lemma is proven. 

We reiterate that there are no other appropriate boundary conditions 
without further assumptions on Tq(p) . Equivalent boundary conditions for 
bounded atmospheres are given in Theorem 2. 

Lemmas 2.1 and 2.2 depend crucially on the entropy assumption (1.24). 
Before proceeding with the proof of Theorem 1, we therefore discuss the 
nature of this assumption. To our knowledge, all thermal profiles 
considered previously in the literature are such that both 


and 


0q -* *# as p -*■ 0 


( 2 . 10 ) 


d0- d0 n 

lim inf (~p Z -rp ) - lim sup (-p Tp ) . 

p->0 ** p-+0 v 


( 2 . 11 ) 


In particular, Staniforth et al . (1985) studied the family 


o (p) - Ap a , A > 0 , (2.12) 

which includes most previously-studied profiles, and they showed that 
linearization validity requires 


a & 2 —k . 


(2.13) 
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Using (2.5) one finds that the potential temperature profile 6 q correspond- 
ing to (2.12) is given by 


^ 0 (p)-«o<p s ) 


while 


( A 1 / -Q+ ( 2 -k ) ~ a+ ( 2-/c ) . 

R a— (2— «) cp P S ; 1 



a > 2 — k 
a - 2— /c 


2 + * _ A 




dp 


R 


(2.14) 


(2.15) 


Clearly both (2.10) and (2.11) are satisfied. 

To complete the discussion of our entropy assumption (1.24), we now 
assert that it is actually less restrictive than properties (2.10), (2.11). 

Claim 2.1 . For thermal profiles satisfying (2.10) and (2.11), one has 

2 +k Q 

lim inf (— p - — ) - + ® . (2.16) 

p-*0 dp 


Hence the entropy assumption (1.24) is satisfied. 

Proof : Suppose (2.16) is false. Then (2.11) implies that 


2+« d0Q 

lim sup (-p ~ ) < <*> . (2.17) 

p->0 P 


Hence there is a constant M > 0 such that 


2 + * d *0 

- P dp" < M for 0 < p < p s 


(2.18) 


Integrating from p to p^, one finds that 
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* 0 <P> < V p s> - i“- 
2 * 



(2.19) 


1 

4 W + p s 

2 K 


where we have used the assumption (1.25b) that k < 1/2. Hence 0 q(p) is 
bounded, contradicting (2.10), and the claim is established. 

So far it has been shown that for any number there is at most one 
square -integrable solution of the VSE: when a solution of the VSE satisfying 
the upper boundary condition (2.8) exists, it is unique . Theorem 1, whose 
proof occupies the remainder of this section, concerns the existence of 
solutions under the additional condition (1.9) at the lower boundary. 

Our theorem is simply an application of a classical result (cf. 
Berkowitz, 1959, Theorem 5.4) in the spectral theory of singular Sturm- 
Liouville problems, namely that the spectrum is discrete below /i — where 

n - lim inf g ^ p) , (2.20) 

c p-*0 4^ (p) 


the function VKp) having been defined in (2.2). In particular, there is a 
countably infinite set of eigenvalues satisfying (1.28) when “ + 00 • To 
prove Theorem 1, we -need only show that defined by (2.20) satisfies 




c 


-H» when T - 0 

CO 

l/^/cRT^) when * 0 . 


(2.21a) 

(2.21b) 


To do so, we prove the following two lemmas regarding the asymptotic 
behavior of o(p) and i>(p) . 
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This is not possible since Tq > 0 for 0 < p < Pg. Hence p dTg/dp has a 
limit . 

Now we prove (2.23) by contradiction. Suppose (2.23) does not hold, 
and define f ® p dT^/dp . Then there are constants S > 0 and p* > 0 such that 


Hence 


i 


e . 





I T 0 <P*) “ T 0 (p)l > s lo S ~ — - 00 as p - 0, 
which contradicts our basic assumption that lim T n (p) is finite 

p-0 

(2.23) must be true and the lemma is proven. 


Lemma 2.4 . 

lim p^(p) - kRT 

p-»0 


Proof : From the previous lemma and the definition (2.2) 

follows that 


where 



g(s) - s o(s) - kRT 


as s -*■ 0 . 


From (2.30) it follows that 


pVKp) " «RT oe (l - 


_D_ 



p 



ds , 


;2.28) 

(2. 29) 

. Therefore 

of iKp) it 

(2.30) 

(2.31) 

(2.32) 


so the lemma is proven if we can show that 
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r -2 

lim p g(s) s ds - 0 . (2.33) 

p-0 J 

P 

This is clearly the case if the integral converges as p ■* 0, so assume it 
does not. Then L' Hospital's rule yields 

r s _ 2 -g(p)P~ 2 

lim p g(s) s ds - lim —z - lim g(p) - 0 , (2.34) 

p-*-0 ' p-»0 — p p-^O 

where (2.31) has been used. This completes the proof. 

The proof of Theorem 1 in case T^ 0 is now immediate, upon rewriting 
(2.20) as 

P 2 <*(P) 

u - lim inf T . (2.35) 

p-»0 4[p^(p)] Z 

Lemmas 2.3 and 2.4 imply that both the numerator and denominator here have 

o 

limits, namely /cRT and 4(/cRT ) , so that the result (2.21b) follows. 

J CO 00 

Finally we treat the case T^ - 0 and complete the proof of Theorem 1 
by verifying (2.21a). To do so we prove the following lemma. 


Lemma 2.5 . For 0 < p < p s> 


RT (p) 

(i) J o(s) ds < — “ 

P 


a( p) _1_ , d l0g 9 _0 v 

(ll) 2 RT n ( P dp } 

P) 0 


0 /, 


Proof: 


Using (1.2b) and (2.2), 


we prove statement 


(i) by 


by parts: 




*-l d# 0 „ 

p dT dp 


-R^p' 4 1 9 q S + (l-«) | s' 1 2 0 Q (s) ds j 


< Rp 


K — 1 


- RT 0 /p , 


where the inequality was obtained since k < 1 and 8q> 0. 
definition (1.2a) of a we then obtain 

<7 P 2 RT 0 , d log 9 _0 . _1_ , d log 9 _0 v 

* 2 > r 2 t 2 p < dp > “ RT q (_P dp > 
thus proving statement (ii). 

From part (ii) of this lemma it is clear that the proof 
is complete if we can show that 


lim inf 

p->0 


RT, 


("P 


d log 9 ( 
dp 


+ 00 


Since Tq -*■ 0 as p -*■ 0, this is certainly true in case 


d log 8 Q 

lim inf (~p 

p-»0 


integrating 


(2.36) 

From the 

, (2.37) 

of Theorem 1 
(2.38) 


dp 


) > 0 . 


(2.39) 
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The only other possibility is that 


lira inf (-p 

p->0 


d log g Q 

dp 


0 , 


since — d log fl^/dp > 0 for 0 < p < p g by static stability. 

d log 6 q d log T Q 
p - . p « . 


But lim (p d log T n /dp) exists by the limit assumption, so 

p-+0 

that in fact 


lim (-p 

p-*0 


d log g Q 
dp 


- 0 . 


Hence, for each e > 0 there exists a p^ > 0 such that 


“ P 


d log g Q 
dp 


< e 


for 0 < p < p + . 


Integrating from p to p^, one obtains 


and therefore 


log 


yp> 

Vp*> 


P 


Mi) 


Now 


Tq < const, x 


K — € 


P 


RT„ ( P 


d log 8 


0 


rt; 


l+K . 

( "P dp } 


(2.40) 

Now 

(2.41) 

(2.40) implies 

(2.42) 

(2.43) 

(2.44) 

(2.45) 


dp 


(2.46) 
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so equation (2.45) implies that 

(2.47) 

(2.48) 

Since k > 1/4 by (1.25a), we may choose e so that 0 < 0. The result (2.38) 
then follows from (2.47) by virtue of the entropy assumption (1.24), thereby 
completing the proof of Theorem 1. 


a o+k 

) > const, x p P (~p ) 


0 - -2(/e - r - e) . 



3 . Bounded atmospheres 

In this section we prove Theorems 2 and 3. The first lemma expresses 
boundedness of the basic state atmosphere as a condition on the static 
stability function cr. 


Lemma 3.1 . Suppose the basic state atmosphere is bounded, i.e., $_ 

E 

defined by (1.18) is finite. Then 


(i) Tq ( p) -*• 0 as p -* 0 

and 

p r s 

(ii) I sct(s) ds - /c$ E - RTg < ® . 

0 


( 3 . 1 ) 


( 3 . 2 ) 


Proof : By the definition of (1.17) of $g(p) and from our basic 

assumption that lim T n (p) exists, it follows that the limit must in fact be 

p-0 u 

zero, for otherwise the integral in (1.17) would diverge at least 
logarithmically . 

To prove part (ii), write (1.17) as 


r s _i 

$ q(p) - VV + R J *0^ s K L ds 


* o (p s> + ; V s) 


5 mm 


s „ d» 0 (s) 
s ~TT~ ds 


]• 


(3.3) 


or, using (1.2b), 

Vp> - W + K ET s - K RT 0 ( P> + i jf "O) ds 


p, 

p 


(3.4) 


From (3.1) one then obtains 
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lim so(s) ds - k lim $ n (p) ~ - RT C < ® , (3.5) 

p->0 ' p^O u ^ s 

P 


completing the proof. 

The following two lemmas are fundamental and will be used extensively. 
The first of these simply states several versions and immediate consequences 
of the variation-of -constants formula for the VSE. 


Lemma 3.2. 
and 0 < p* s p g . 


Suppose that X(p) satisfies the VSE (1.8). Let 0 < p < p 

J 

Then 


(i) 


where 


X(p) - X(p*) + a<Kp) - n 



P 


[iKp)-^(s)] X(s) ds 



a(s) ds 


(3.6) 


(3.7) 


1 

o 


dp 


+ n X(s) ds - const. 


(3.8) 


Here ij> - ^(p) depends only on p* ( while a depends on p* but not on p. 
(ii) If X satisfies the LBC (1.9) then X(pg) * 0 and 


X(p) 



X(P) 3 X(p) / X(p s ) 

? s 

V>(p) - A» J [iKp) - V>(s)] X(s) ds 

P 


(3.9) 


(3.10) 


satisfies 


0 as p -* 0 then 


(iii) 


If if 
a dp 


n r\ 




1 dX 
o dp 


- /i J X(s) ds 
0 


(3.11) 


and 


X(p) 


P* 


- X(p*) + MV>( P) j X(s) ds - n j [*(p) - rf>(s)] X(s) ds (3.12a) 


P* s 

X(p*) + M I [1 X(p') dp ' J ct(s) ds . 

p 0 


(3.12b) 


(iv) Ifif->0asp->0 and if X satisfies the LBC (1.9), then 
a dp 


_P< 

H! 


f 

J 


X(s) ds , 


(3.13) 


where X is defined in (3.9). 


Proof : By differentiating (3.6), dividing by o(p), and 

differentiating the result, one obtains the VSE. Thus (3.6) defines a 
solution of the VSE. The same argument can be applied in reverse by virtue 
of the continuity assumptions, so that the VSE is in fact equivalent to 
(3.6) and part (i) is proven. 

In part (ii), if X(p g ) - 0 then the LBC would imply dX(p g )/dp - 0, 
hence the VSE would imply X(p) ■ 0. To verify (3.10), evaluate (3.8) at 
p = p g and use the LBC to find 



X(p s ) - ^ 


r 


p* 


X(S) 


ds . 


a - 


(3.14) 
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Now substitute (3.14) into (3.6), choose - p s , and divide by X(p g ) to 
obtain (3.10). 

To verify (3.11) in part (iii), let p ■+ 0 in (3.8), so that 

P* 

a - j X(s) ds , (3.15) 

0 


then substitute (3.15) into (3.8) to obtain (3.11). Equation (3.12a) 
follows by substituting (3.15) into (3.6). Equation (3.12b) follows most 
easily by multiplying (3.11) by o(p) and integrating the result from p to 
P*> 

Part (iv) is verified by equating (3.14) and (3.15), and then dividing 
by X(p g ) . 


Lemma 3.3 . Suppose that the basic state atmosphere is bounded, and 
that X satisfies the VSE (1.8) and LBC (1.9). Then 

(i) X has finitely many zeros. 

(ii) If fi < 0 then in fact X has no zeros. 

(iii) Both X and ^ ^ have limits as p •* 0. 


Proof : The zeros of any solution of a singular Sturm-Liouville 
equation can accumulate only at the singular point. First we show that if X 
has any zeros at all, then fj. > 0, thereby proving statement (ii). Denote by 
p^ the location of the first zero to the left of p<,, so X defined in (3.9) 
is such that X(p^) - 0 and X(p) > 0 for p^ < p < p g . Setting p - p^ in 
(3.10) then gives 


1 + 


RT„ 


V>(P 


r s 

L ) - M J ItKPj.) -*(»)] X(s) da , (3.16) 
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from which it follows that n > 0 because V’(p^) > 0 and the integrand is 
positive . 

To prove (i) , suppose n > 0. Let p^ and p^ be two successive zeros of 
X(p) , with 0 < < Pj < Pg. Choosing p* - Pj in (3.6) gives 


where 


X(p) - a ^(p) 



ds , 


(3.17) 



(3.18) 


Now either X(p) > 0 or X(p) < 0 for p^ < p < Pj . In the former case, since 
H > 0 we find from (3.17) that 


0 < X(p) < a ^(p) for p t < p < pj , 
hence a > 0, while in the latter, 


0 < - X(p) < - a V>(p) for p t < p < pj , 
and a < 0. Combining, we have 


(3.19) 


(3.20) 


0 < | X(p) | < | a | ^(p) for < p < Pj , (3.21) 

and a * 0. 

Setting p - p^ in (3.17) gives 


0 


P f J 

a “ M J llKPj.) _ V»(s)] X(s) ds , 


(3.22) 
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so that 

P (j 

|a|^(p.) - n J Wp.) - V»(s) ] |X(s) | ds 

Pi 

h 

< M J [V’(P i ) - iK S)] | a| ^(s) ds 

P i 

fj 

< M J ^( ?i ) | a | tf(s) ds (3.23) 

p i 


where (3.21) has been used. Therefore 




Pi 


(s-p^ a(s) ds 


< M 


V 

I sa(s) 

Pi 


ds 


(3.24) 


Summing over all pairs of successive zeros of X, it follows that the total 
number n of zeros for 0 < p < Pg satisfies 
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n-1 < JU I sa(s) ds . (3.25) 

0 

The integral is finite by Lemma 3.1, hence part (i) is proven. 

Finally we prove part (iii). Our continuity assumptions on a 
guarantee that every solution X of the VSE is continuous for 0 < p ^ Pg. 

fP* 

Since X also has finitely many zeros it follows that lim I X(s) ds exists 

p->0 •'p 

for every p^ > 0. Therefore (3.8) implies that ^ has a limit as p ■+ 0. 
Now choose p^ such that X(p) does not change sign for 0 < p < p^. Then it 
follows from (3.8) that dX/dp has at most one zero for 0 < p < p 
Therefore X is monotone in some neighborhood of p - 0, and hence X has a 
limit as p -*• 0. This completes the proof of Lemma 3.3. 

The most difficult part of the proof of Theorem 2 is showing that 
statement (iii) of the theorem implies statement (ii) of the theorem. We 
therefore begin by proving a lemma which gives several consequences of 
statement (iii). 


Lemma 3.4 . Suppose the basic state atmosphere is bounded, and that X 
satisfies the VSE (1.8), the LBC (1.9), and J ^ - 0 as p - 0. Then 


(i) 

M >0 

(3. 

.26) 

(ii) 

pX(p) -*• 0 as p -» 0 

(3. 

• 27) 

(iii) 

If X has any zeros , then 





where p^ denotes the first zero of X(p) to the right of p - 0. 
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(iv) Given any e > 0, there is a p Q > 0 such that 

(i j s ct(s) ds < £ if 0 < p < Pq . (3.29) 

0 

If € < 1 and 0 < p < Pg, then 

P 

j |X(s) | ds < ^ p|X(p)| . (3.30) 

0 

Proof : To prove part (i), suppose that n i 0. We then have 

X(p) ■ X(p)/X(p s ) > 0 for 0 < p ^ p g by Lemma 3 . 3 ( ii) , hence the right-hand 
side of (3.13) is nonpositive. This is impossible since the left-hand side 
is positive. Therefore n > 0. 

To prove part (ii), observe that since X has finitely many zeros 
by Lemma 3.3(i), there is a p* > 0 such that X(p) * 0 for 0 < p < p*. Since 
fi > 0, for this choice of p^ in (3.12a) we find that for 0 < p < p^, 
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where we have used (3.8) and Lemma 2.5(i). Therefore 


0 < p |X(p) | < p|X(p*)| + 1 a | RTq(p) 

■+ 0 as p -* 0 , 

(3.33) 

by (3.1), and part (ii) of the lemma is 

established. 


Now let p^ be as defined in statement (iii). Then setting 

P* " Pi 

(3.31) yields 




?1 

?i 



|X(p) I < n <7(s) ds 

|X(s)| ds , 

(3.34) 

P 

0 



and integrating (3.34) from zero to p^ leads 

to 


f 1 r f 1 i 

1 < M «^(s) ds 

) L J 

dp 



0 p 




/ p r 1 

P 1 

p f l l 


- M \ P o(s) ds 

1 + 

p o(p) dpi- 


L J 

P 

0 

J J 

0 


- M P <*(p) dp - ^ 

lim p 

> <r(s) ds , 

(3.35) 

J 

0 

p-+0 

J 

P 



provided the indicated limit exists. In fact the limit is zero since, by 
Lemma 2.5(i) , 


Pi Pg 

0 < p J cr( s) ds < p J cr(s) ds < RTq(p) 0 as p -► 0 . (3.36) 


This proves statement (iii) . 
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The first part of statement (iv) follows immediately 
Lemma 3.1(ii). If e < 1 then, in particular. 


fo 

s ff(s) ds < 1 . 

0 


(3.37) 


If X has any zeros, then (3.28) and (3.37) together imply that pg 
Hence X(p) * 0 for 0 < p < pg. But from (3.11) we have 


dp'] , (3.38) 

and integrating from zero to p and using (3.27) gives 


% <px> - * + <p»xi §> 


X + <p»)U< j X(p') 
0 


p p s 

pX(p) - J X(s) ds - n J s a( s) £ J X(p') dp'] ds . (3.39) 

0 0 0 


Since X(p) does not change sign for 0 < p < p^ we therefore have on 
interval 


p p s 

p|X(p) | - J |X(s) | ds - n J sa(s) [ j |X(p')| dp'] ds 
0 0 0 


> 


I 


|X(s)| ds 



sa(s) 





0 


|X(s)| 


ds 


from 


Pr 


this 
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> (l-o | |X(s) | ds , (3.40) 

0 

where the last inequality follows from (3.29). This completes the proof of 
Lemma 3.4. 

We are now ready to prove Theorem 2. Referring to the statements in 
the theorem numbered (i) through (v) , we shall prove the theorem by showing 
that 

(iii) -> (ii) -> (i) -> (iii) 

(i) <-> (iv) 

(i) <-> (v) 

Proof that (iii) -> (ii) : Let e and Pq be as stated in Lemma 4(iv) , 

and let 0 < p < p* < Pq. From (3.12b) and (3.30) it then follows that 


|X(p)| < |X(p*)| 




ct(s) 


|X(s) | ds . 


(3.41) 


By virtue of the Gronwall inequality (cf. Coddington and Levinson, 1955, 
p. 37, or Hale, 1969, pp. 36-37), it follows from (3.41) that in fact 


|X(p)| < |X(p*)| jl + ^ j s a( s) exp J ^7 J P' CT (P') d P'] ds } • (3.42) 


The integrals here may be bounded in terms of e, using (3.29), as follows: 


I X ( p ) | < |X(p*)| 



-JL. 

1—e 



P 


s a(s) 


exp( I^7 ) 
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< |X(p*)| {l + expC^)} < 00 . (3.43) 

But X(p) has a limit as p ■* 0 by Lemma 3.3(iii), hence (3.43) implies that 
the limit is finite. 

Proof that (ii) -> (i) : Every solution of the VSE is continuous, and 
therefore bounded, for p^ < p £ p^ with p^ > 0 arbitrary. Statement (ii) 
implies that in fact this is true for 0 ^ p < p g . Therefore X is square- 
integrable . 

Proof that (i) -> (iji) : Lemma 3.3(iii) guarantees that 

w ■ u “ Si (3 - 44) 

p-*0 ^ 

exists. We will show that (i) -> (iii) by proving the contrapositive. 
Suppose that W ** 0, so that there is a p + > 0 such that 

ij S 1 > 2 w > 0 for 0 < p < p* • < 3 - 45 > 

From Lemma 2.1 it then follows that 

|^| > \ W c iP - 3/2 for 0 < p < p* . (3.46) 

Upon integrating from p to p^, this yields 

| X(p) - X(p*)| > W Cl (p _1/2 - p* 1/2 ) for 0 < p < p* . (3.47) 


It follows that 
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r s 2 r 2 

X z (s) ds > X^(s) ds 


l* 


> j |X(s)— X(p*)| 2 ds - [ X 2 (p*) ds 


J* 

J “ '***; 

P 


> W 2 c 2 | (s 1/2 - p* 1/2 ) 2 ds - (p*- p) X 2 (p*) 

P 

-►oo as p -*■ 0 . (3.48) 


Proof that (i) -> Civ) : We have already shown that (i) implies both 

(ii) and (iii). But (ii) and (iii) together imply (iv). 

Proof that ( ivl -> ( i") : Statement (iv) implies that lim X(p) - 0 or 

p-»0 

lim — ^ - 0, since Lemma 3.3(iii) states that both of these limits exist 

p-0 a dp 

for every solution of the VSE. In the first case, (i) follows trivially. 
In the latter case (i) also follows since we already showed that 

(iii) -> (i). 

Proof that ('ll -> (v) : This is trivial since (i) -> (ii) -> (v) . 

Proof that (v) -> (i) : We will show that (v) -> (iii), by proving the 
contrapositive. As shown previously, the negation of (iii) implies (3.47). 
Multiplying (3.47) by p^ gives 

|p 5 X(p) - P 5 X(p*)| > Wc^p 5-172 - P 5 P ; 1/2 ] , (3.49) 

and therefore, letting p -*■ 0, we have 

lim inf p 5 |X(p)| > Wc t > 0 , 

p-0 


(3.50) 



for 0 < S < 1/2. This completes the proof of Theorem 2. 


/, n 


'-r 


Proof of Theorem 3 : Part (i) follows directly from Theorem 1 and 

Lemma 3.1(1). We do not prove part (iii), but only point out that for 
limit-point Sturra-Liouville problems with totally discrete spectra, the 
facts about completeness, orthogonality, and number of zeros of the 
eigenfunctions are still the same as for regular Sturm- Liouvi lie problems; 
cf. Coddington and Levinson (1955, Chapter 9, Theorem 3.1 and Problem 1). 

We now prove part (ii). Positivity of the eigenvalues was already 
stated and proven as Lemma 3.4(i). To obtain the eigenvalue bounds we 
use the fact that the eigenfunction X n corresponding to eigenvalue n n has 
precisely n zeros. In particular, Xg(p) ■ Xg(p)/Xg(p s ) > 0 for 0 < p ^ p g . 
Since Mq > 0. it then follows from (3.11) that dXQ/dp < 0 for 0 < p ^ p g , so 
in fact Xq(p) > X Q (p g ) - 1 for 0 < p < p g . From (3.13) we then have 

> VS ’ <3 ' 51) 

proving the upper bound in (1.31a). 

To obtain the lower bound, use (3.13) and Lemma 3.4(ii) to write 



(3.52) 
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Now from (3.11) and (3.13) we have for 0 < p < p g 


dX, 


1 r_o 

a dp J '“'0 

0 


J X n (s) ds 


P S 

< /* 0 j V s ) ds 


RT S * 


(3.53) 


Introducing (3.53) into (3.52) and using Lemma 3.1(11), we get 


1 ? S 

rt: < “o ♦ rt: J s ° w 


^ £ 

**0 RTg ’ 


(3.54) 


which completes the proof of (1.31a). 

The bound (1.31b) for n 2 1 is actually a refinement of (3.25). If we 

denote the first zero of X by p. and the last by p . we have from (3.24) 

n 1 n 

that in fact 




n— 1 < fi I s a(s) ds . 


(3.55) 


Adding inequalities (3.28) and (3.55) then gives 



42 


r 

n< ^n J 


s <j(s) ds 




< H s o(s) ds 


(3.56) 


and (1.31b) follows from Lemma 3.1(ii). 

Inequalities (1.32a,b) follow directly from (1.31a,b). We remark 
that, while the bounds in (1.32a) usually provide a good estimate of the 
barotropic equivalent depth as pointed out in the Introduction, the bounds 
in (1.32b) become progressively worse as n increases. This is clear from 
the derivation of (3.24). Better bounds can be obtained if one has some 
information about Tq(p) . 

It remains to prove statement (iv). If p 1 denotes the first zero of 
X n , it follows from (3.11) that X fl (p) is monotone for 0 < p < p^, hence 

|X n (p)| > |X n (p*)| > 0 (3.57) 


for 0 < p ^ p^ < p^. Again using (3.11) we then have 


I dX n (p) 

I dp 


P 

/* n <Kp) \ l x n ( s >l ds 
0 


> 4 n |X n (p*) | pa(p) 


- 1/2 


>M n |X n (p*)| C] p 


(3.58) 
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for 0 < p < p* < p^, where we have used Lenuna 2.1. This completes the 
proof. 
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4 . Examples 

In this section we give four examples that illustrate some of our 
theoretical results. In the first two examples, we specify the asymptotic 
form of Tq(p) near p - 0 and we apply the Frobenius (power series) method 
(e.g., Coddington and Levinson, 1955, Section 4.8) to determine the 
asymptotic behavior of all solutions of the VSE. These asymptotic results 
depend neither on Tq(p) away from p - 0 nor on the lower boundary condition. 
In Example 1 the atmosphere is bounded and the asymptotic results allow one 
to verify directly the equivalent boundary conditions of Theorem 2. In 
Example 2 we have lim T n (p) * 0 and we show that the spectrum is continuous 

rwf> u 
r - 

for ^ a result which complements Theorem 1. 

In the third and fourth examples we specify T Q (p) throughout the 
domain, and we obtain closed-form solutions of the VSE satisfying the lower 
boundary condition. In both examples Tq(p) -* 0 as p •+ 0, but slowly enough 
that the atmosphere is unbounded. These two examples are used as counter- 
examples showing that Theorem 2 is not true when its stated hypothesis of a 
bounded atmosphere is replaced by the somewhat weaker hypothesis that 

Tq(p) -► 0 as p -► 0. In particular, Example 3 shows that it is possible for 

1 dX 

non- square -integrable solutions of the VSE to have — — -► 0 as p 0, while 
Example 4 shows that it is possible for square- integrable solutions (eigen- 
functions) to be unbounded. 

Example 1 . Suppose that near p - 0, T Q (p) is an analytic function 

of p , 

CO 

T n (p) - S r p €J , (4.1) 

u j-o J 

where we assume that 

0 < e < ac , - 0 , r l > ^‘ 


(4.2) 
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It follows that our basic assumptions are all satisfied near p - 0, and that 


Ta, - lim T 0 (P) - 0 . (4.3) 

p-*0 


Hence the spectrum is totally discrete. 

From (4.1) and the definition (1.2c) of a, it follows that near p - 0, 
a has the form 

00 

a(p) - Rp € ~ 2 S [* - «<j+l)] r . . p €j . (4.4) 

j-0 J+1 

By taking - 0 for j > 2, this family reduces to the family (1.14,1.15) 
considered by Staniforth et al. (1985), not including the limiting cases 
a - 2 — k and a - 2. 

Using (4.4) we write 

where 

* 1 

B(q) - S b.q J ; (4.5b) 

j-0 J 

the coefficients b^ can be calculated from the coefficients in (4.4). Now 
substituting (4.5a) into the VSE (1.8), and making the coordinate transfor- 
mation 

q - P £ , (4.6) 

the VSE becomes 

2 

q 2 + qC(q) jg + D(q)X - 0 , (4.7a) 

dq " 

where 


C(q) - q ^ (log q 1/£ B(q)) 


(4.7b) 
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- ^2 • (4 - 7c) 

The Frobenius method applies since the coefficients C and D are 
analytic functions of q. Therefore the solutions of the VSE have the form 

X(q) - q^[l + 2^ aj q j ] (4.8) 

in a neighborhood of p - 0, where the series converge uniformly and can be 

differentiated term-by-term. The exponent f) and coefficients a^ are obtained 

by substituting (4.8) into (4.7a), using the power-series expansions for 

C(q) and D(q), and then equating coefficients of like powers of q. Using 

(4.6) one finds that the regular solution X^(p) - X^(p;/i) and the 

( 2 ) ( 2 ) 

singular solution X (p) - X (p;/i) are given near p - 0 by 


X (1) (p) - 1 - V p£ + S a i 1) P <J • 
U £ j-2 J 


(4.9a) 


X (2) (p) - p £ 1 [l + S a< 2) p ej l 

L j.i J J 


(4.9b) 


where the coefficients and a j^^ can be expressed in terms of the r ^ in 

(4.1). 

It is a simple matter to verify that X^^ satisfies each of the 

( 2 ) 

boundary conditions of Theorem 2, whereas X v satisfies none of them. In 
particular the n^ eigenfuncion is given by 


x n (p) - X (1) (p;/i n ) , 


H being the n C ^ 


eigenvalue . 


(4.10) 
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Example 2 . Suppose that near p - 0, T n (p) is an analytic function 


of p, 


T 0 (P) “ 2 r p J , r Q > 0 . 

j-0 J 


(4.11) 


If rj - 0 for j > 1, this reduces to the familiar atmosphere with an 
isothermal top. Our basic assumptions are all satisfied near p - 0, and 


T co " lim T 0 ( P> " T n * 0 

p->0 


(4.12) 


Th i the atmosphere is unbounded and from Theorem 1 it follows that the 
por .on of the spectrum for 




is discrete. We will show that the spectrum is continuous for p > n 
From (4.1) and the definition (1.2c) of a, we have 


CO 

o( p) “ Rp~ 2 2 (/c-j) T pJ 


refore 


j-o 


7vT - p 2b < p3 ■ 


(4.14) 


(4.15a) 


where B is the analytic function 


B(p) - S b p J 

j-0 J 


(4.15b) 


whose coefficients can be expressed in terms of the t ^ . 
written as 


The VSE can then be 


P 2 + P C(p) g + D(p)X - 0 , 
dp 


(4.16a) 
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where 


and 


C(p) = 


1 dp 2 B 
pB dp 



(4.16b) 


(4.16c) 


are analytic. Application of the Frobenius method results in the following 
three cases. 

If fi < fi we find that near p - 0 the two linearly independent 
solutions are of the form 

& 00 

X (1) (p) - X a) (p;/i) - p + [1 + 2 a. (1) pj] (4.17a) 

j-1 J 

X (2) (p) - X (2) (p;,i) - p^"[l + 2 a (2) p j ] (4.17b) 

j-1 J 

where 



(4.17c) 


Clearly X^^ is square- integrable and X^ 2 ^ is not. For any eigenvalue 
< M c it follows that the corresponding eigenfunction is 


X n (p) - X (1) (p; Mn ) . 


(4.18) 


The actual number of eigenvalues /i < m c> if any, depends upon the thermal 
structure all the way to p - p<,, as well as upon the lower boundary 
condition. This number is always finite. 

If/* — a* , we find instead that 


x ( 1 ) (p) 


p- 1/2 ll 


+ 2 
j=l 




(4.19a) 
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X (2) (p) - P _1/2 (log p) [1 + S aJ 2 Vl . (4.19b) 

j“l J 

It is clear that is not square - integrable . Since |X^(p)| > |X^(p)| 

( 2 ) 

for all p sufficiently small, it follows that X 7 is not square- integrable 
either, nor is any linear combination of X^ and X^ 2 ^ . 

If n > n c , two linearly independent solutions are 

CO 

X ( 1 ) (p) - p - 1 / 2 [cos (7 log p)] [1 + <7 a. ( 1 ) p j l (4.20a) 

j -1 J 

and 

X (2) (p) - p~ 1/2 (sin (7 log p)] [1 + a a. (2) p J ] (4.20b) 

j -1 J 

where 

7 - 2 J”j{“ ” 1 • (4.20c) 

Observe that X^(p) is not square- integrable, for 

J S 1 cos 2 (7 log s) ds - — 2 log p - 2 ^” sin (27 log p) , (4.21) 

P 

( 2 ) 

which is unbounded as p -» 0. Similarly it follows that X v is not square - 
integrable, nor is ,any linear combination of X^ and X^ 2 \ 

Thus there are no square- integrable solutions of the VSE for n > n . 
This means that the spectrum is continuous for n > for every such p 

there is a unique solution of the VSE satisfying the lower boundary 
condition (1.9), but the solution is not square - integrable . However, one 
can still expand arbitrary square - integrable functions in terms of the 
totality of these solutions for n > n c and the eigenfunctions (if any) 
for (i < in a manner analogous to the Fourier integral; cf. Coddington 
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and Levinson (1955, Chapter 9, Theorem 3.1), and Staniforth et al. (1985, 
p. 346). Thus if one is interested in spectral expansions, then no upper 
boundary condition should be imposed when p > We reiterate that 

continuous spectra cannot arise for bounded atmospheres. The continuous 
spectrum in this example is spurious in the sense that the primitive 
equations on which the VSE is based do not apply to unbounded atmospheres, 
due to the traditional shallowness approximations. 

Since the atmosphere is unbounded in this example, Theorem 2 does not 
apply. In particular, for all p > 0 it follows from our three sets of 
asymptotic formulas that both and X^^ satisfy 

J § - 0 as p -* 0 . (4.22) 


Thus the boundary condition (4.22) is generally not well posed, because 

it is not strong enough to eliminate the non- square- integrable solutions 

when u < u . 

^ 'c 

The boundary condition (4.22) is important because it is derived 
from the commonly invoked boundary condition (1.6) that w 0 as p 0 in 
the primitive equations. Theorem 2 implies that (4.22) is well-posed for 
bounded atmospheres, while Example 2 shows that it is generally not 
well -posed when T^ * 0. In the next example T^ - 0 but the atmosphere is 
unbounded. We will, see that the boundary condition (4.22) is still not 
well-posed. 

Example 3 . Let > 0 be arbitrary, and let e > 0 be an arbitrary 
parameter. Consider the VSE problem for 


RT 

tf(p) - — | (1 - « log 


JL- x-2 


2 

y s p s 

RT„ 


- € log 


( 


) • (4.23) 
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f 

i 
| 

From (1.2b) it follows that the temperature profile Tg(p) corresponding to 
any a(p) is given by 

T 0 (P) " P* TgPg* + R -1 | s 1_K o( s) dsj . (4.24) 

P 

t 

Upon substituting (4.23) into (4.24) it is easily verified that our basic 
assumptions in Sec. 1 are all satisfied. 

It can also be established that 


T - lira T A (p) - 0 , (4.25) 

p-0 ° 

as follows. For a given by (4.23) it is clear that the integral in (4.24) 
diverges as p -*• 0. Therefore L' Hospital's rule applies to give asymptoti- 
cally 

vp> - i [ ) - p - o 

v — K p ' 

1 2 , . 

- p a(p) 

- T_(-k£ log ) 1 , (4.26) 

b P S 

from which (4.25) follows immediately. By Theorem 1 the spectrum of the VSE 
is therefore totally discrete. Nevertheless the basic state atmosphere is 
unbounded because from (4.26) we have that 

l 

p 


RT Q (s) 


r ds 

Ke ' — s log ~ 

° TV 
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RT 

- 7~ log(-log -J- ) , (4.27) 

P s 


so from (1.18) it follows that $q(p) H ® as p 0. 

By differentiation it can be verified that the function 


where 


X(p) 



-i//e 

(1 - e log ) 

P S 


u 


satisfies the LBC and VSE with 


_ !sZa 

RT C 


A* - 



(4.28) 


(4.29) 


(4.30) 


On the other hand, this X is not an eigenfunction because it is clearly not 
square- integrable . However, one finds that 


idU'-u-' log -j- r 0, 0/ ‘ • (4-3D 

o 

Therefore (4.22) holds whenever the parameters are such that e < v. Once 
again we see that the boundary condition (4.22) is too weak to exclude 
solutions of the VSE which are not square- integrable. This example also 
shows that Theorem 2 does not hold under the hypothesis that Tq -* 0 as 
p -*• 0, which was shown by Lemma 2.1 to be weaker than the stated hypothesis 
that the atmosphere be bounded. 

Notice also that if e > u, this example shows that it is possible for 
a function to satisfy the LBC and the VSE with \x < 0. 

We remark that the parameters a <, and e in this example can be 
chosen in such a way that T Q (p) is arbitrarily close to being constant on 
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any interval < p < pg with p^ > 0. To see this, choose a g so that 


dT Q /dp - 0 at p - p s , i • e . , 


CT s ’ * R V S 2 


(4.32) 


Then (4.23) reads 


-2 


<r(p) - kRT-p (1 - « log ~T~ ) (1 ~ T log ) 


-B— n ~ 2 n - 4 


-2 


- /cRT s p ‘[1 + 0(0] 


(4.33) 


Therefore a(p) can be made arbitrarily close to the isothermal static 

-2 

stability function /cRTgP over any interval away from p - 0, by choosing 
e sufficiently small. For such a choice of parameters, it follows that 
T 0 (p) - Tg from the surface up to a point p^ arbitrarily close to p - 0, 
after which (4.26) applies and T^(p) drops to zero. 

It is well known, and also follows from Example 2, that the VSE for 
isothermal atmospheres has a continuous spectrum. Example 3 shows that by 
modifying the isothermal atmosphere in a neighborhood of p — 0, a totally 
discrete spectrum can be obtained. 

This example was created by choosing ^ ^ (judiciously) to be of the 
form (4.31), then differentiating to obtain -/iX, then differentiating once 
more to obtain ~n dX/dp. The LBC then determines n, and a is recovered upon 
dividing dX/dp by ^ • The following example is constructed in the same 

1 V 

way, but from a different choice of — — . As in Example 3, we will have 
T — 0 but an unbounded atmosphere. The example will exhibit an 

oo 

eigenfunction which is unbounded. This is in contrast with Theorem 2, which 
states that for bounded atmospheres all eigenfunctions are bounded. 

Example 4 . Let a g > 0 be arbitrary, and let a parameter S satisfy 


0 < 8 < min( 


!s£s 

RT„ 


. 1 ) • 


(4.34) 



54 


Define for given ct c and S the positive parameter e , 

J 

€ 3 - 5) • (4 - 35) 


and consider the VSE problem for 


RT, 


^(P) 


(1 - 


• 108 p s ’ 


-2 


r Vs 

RT„ 


( ^ - s > 108 1 : ] 


(4.36) 


As in Example 3 it can be verified that the basic assumptions are satisfied, 
that Tq(p) 0 as p ■* 0, and that the atmosphere is unbounded. In particular 
the spectrum is totally discrete. One can make the temperature profile 
corresponding to (4.36) nearly isothermal as in Example 3, by the choice 
(4.32) of a„ and by choosing S sufficiently close to k. For the choice 

k) 

S - 1/2, (4.36) is identical to the static stability function (4.23) of 

2 

1 

Example 3 with e - — 2 • 

S 

By differentiation it can be verified that the function 


X(p) ■ (1 - c log ~T~ ) 

P S 


($“«)/< |“ 5 _ 


< - 5) 108 i; 


(4.37) 


satisfies the LBC and VSE with 


Clearly we have asymptotically 


l-S 

<“RT S 


(4.38) 


X <P> - < < ^ " S) (1 " ‘ 108 )SA 


( 4 . 39 ) 
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so X(p) ■+ ® as p -* 0. However, this X is square- integrable because 


| (- log x) 25/e dx - + 1) , (4.40) 

0 

where T is the gamma- function , which is finite. Thus X is an unbounded 
eigenfunction. In fact, since X has no zeros, p given by (4.38) is the 
lowest eigenvalue. 
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5 . Concluding remarks 

In this article we have presented an analysis of the vertical 
structure equation for quite arbitrary thermal profiles, using classical 
techniques for singular S turm-Liouville problems. We have argued that the 
vertical structure equation is assuredly meaningful only for bounded 
atmospheres because the primitive equations from which it is derived apply 
only to bounded, in fact shallow atmospheres, owing to the traditional 
shallowness approximations. For bounded atmospheres we have shown that the 
spectrum of the vertical structure equation is always totally discrete, that 
the barotropic equivalent depth corresponding to the lowest eigenvalue 
satisfies the inequality RTg/g < Hq < Jcfc^/g independently of the thermal 
profile, that the eigenfunctions are all bounded but have unbounded first 
derivatives, and that any one of a number of upper boundary conditions, 
including conditions which follow from requiring mass or energy conservation 
in the linearized primitive equations, are equivalent and well posed. These 
results are in stark contrast to the situation for unbounded atmospheres, 
for which we have shown that typically the spectrum is partly continuous and 
the eigenfunctions, if any, are square -integrable but unbounded. 

The primitive equations, with the traditional shallowness approxima- 
tions, are only a model of the atmosphere. Our results pertain to this 
model and not necessarily to any real atmosphere. Under the model 
assumptions there is a vertical structure equation and its spectrum is 
totally discrete. A primitive equation system without shallowness 
approximations would be a more realistic model of the atmosphere. It will 
be interesting to determine whether such a model has a vertical structure 
equation and, if so, the extent to which results on the spectrum of and 
boundary conditions for the classical vertical structure equation carry 
over. These issues will be studied in a forthcoming paper. Work on 
numerical treatment of the vertical structure equation and the primitive 
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equations themselves, in light of the equivalent boundary conditions of 
Theorem 2, is also underway. 
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